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1 Introduction 

Planar = 4 super Yang-Mills is integrable. This means that this theory will be solved 
exactly. Recent advances on two-loop scattering amplitudes [1-4] lend renowed credence 
to the long standing hope, that compact analytic formulas describing nontrivial multiloop 
scattering processes can be found [5]. 

Compared to the now relatively mature problem of the spectrum of anomalous dimension 
of local operators, the problem of computing the S-matrix is still rather in its infancy. In 
the spectrum problem, integrability was the key to unlocking the theory [6, 7]. Integr ability 
made it possible to calculate, to any desired accuracy, quantities with a nontrivial dependence 
on the t'Hooft coupling A = ^ym-^c arbitrary finite value of the coupling [8-11]. On the 
other hand, the recent advances in S-matrix computations did not really exploit integrability, 
beyond the duality with Wilson loops which was used to afford some simplifications. The 
methods used, ranging from physical intuition regarding collinear limits and Regge limits, to 
an important innovation in the handling of transcendental functions called the "symbol" , are 
not specific to integrable theories. 

The fact that integrability was not central should be viewed as a good thing, as this 
suggests a potential for leaking these developments beyond the realm of A/" = 4. On the 
other hand, it appears likely that a good grasp of the symmetries will become increasingly 
necessary as we progress toward the depths of this theory. Contact with integrability in 
scattering amplitudes was made from the string side at strong coupling [12]. Arguably, it 
still remains elusive at weak coupling, despite at best some glimpses at it through the OPE 
approach to computing Wilson loops [13-15]. There may also some constrains, most notably 
Regge limits [16, 17] or self-crossing relations [18], which although not obviously related to 
the Yangian symmetry, maybe have not yet been fully exploited. 

Yangian symmetry at loop level appears subtle. Most generators of the Yangian are 
known to require modification at loop level, including the ordinary superconformal and dual 
superconformal symmetry. They acquire some right-hand side when acting on amplitudes: 



Here x and 9 are the so-called dual coordinates, explained below. This equation is not an 
"anomaly," in the sense that the right-hand side 5An can be absorbed into an appropriate 
modification of the symmetry generators on left-hand side, without changing the algebra. 
The deformation is well-understood at 1-loop [19-22], but in principle a new formula has to 
be re-derived at every new loop order. We stress that the problem is deeper than infrared 
divergences. The origin of the problem, as understood in these works, lies in the notion of an 
asymptotic particle state prepared at infinity not being conformally invariant. A conformal 
transformation will cause the external state to radiate some collinear particles. This nontrivial 
affects even canonical infrared finite quantities such as the ratio of NMHV to MHV amplitudes 
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[19]. In fact, as stressed in [23], the whole fact that the MHV remainder function in A/" = 4 
is not constant is such an effect. 

In this paper we propose enlarging scattering amplitudes into a bigger object. The idea 
is that within this bigger class of object, no right-hand sides would appear and the naive form 
of the generators should be the correct one, at least for some part of the Yangian algebra. 
Specifically, we propose that a generalization of scattering amplitudes exists, that depends on 
4n new Grassmann variables Xi, which can be viewed as constrained 9's, such that the dual 
superconformal symmetries become 



for any infrared finite observable, such as the MHV remainder function or the ratio function. 

The motivation for this new object comes from the Wilson loop side of the scattering- 
amplitude/ Wilson loop duality. MHV scattering amplitudes are dual to bosonic null polygonal 
Wilson loops [25-28] , and other helicity amplitudes are dual to supersymmetric null polygon 
Wilson loops [29-33]. While the duality requires only a Wilson loop in a chiral half {xi,9i) 
of the superspace, the proposed generalization is a Wilson loop in the full {xi,6i,9i) M = A 
superspace. The previously missing half of the variables, restores the chiral half of the dual 
supersymmetries which was broken.^. 

As a main application, we will see that the proposal leads, when combined with explicit 
component expressions for the super- Wilson loop, to a concrete and remarkably powerful 
method (by today's standards; hopefully, less remarkable by tomorrow's standards) for com- 
puting any derivative of a multiloop scattering amplitude. The main result of this paper is 
an integral representation for the total differential of the two-loop MHV n-point amplitude, 
and its symbol; for n > 7 this was not previously known. 

This paper is organized as follows. In section 2 we review the basic notation and formalism 
we will be using. This includes a six-dimensional notation recently developed by Weinberg, 
which can be useful in an arbitrary conformal field theory, as well as a brief review of symbols. 
In section 3 we state a precise conjecture regarding the supersymmetric Wilson loops, and 
work out its immediate implications for MHV amplitudes. In section 4 we apply it to a general 
two-loop n-point MHV amplitudes and explain how it leads to a formula for the derivative 
of this object. Section 5 describes several consistency checks which we have applied to this 
result. We finish with some conclusions. The symbol of the n-point amplitude is reproduced 
in Appendix A, and in two other appendices we describe a general method for computing the 
symbol of a rational integral as well as the reduction to two-dimensional kinematics. 

^There seems to be no obvious connection with the recent non-chiral formulation of scattering amphtudes 





[24]. 
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2 Notations and preliminaries 



2.1 Six-dimensional notation 

In this paper we will be will be concerned with finite quantities in a conformal field theory, 
and expect conformal invariant answers. In this section we would like to describe a formalism 
which allows to maintain conformal invariance manifest at every stage, in any conformal field 
theory, following a recent suggestion by Weinberg [34]. We found it quite useful in practice. 

We will be working in coordinate space. The usual spacetime coordinates , where 
fi = 1 . . .4, single out a preferred Poincare subgroup of the SO(4,2) conformal group. To 
avoid singling out such a subgroup, it is advantageous to package the coordinates into a null 
6- vector 

^ x^" 

X+ \= 11 (2.1) 



X' = 



X- 




where X's are identified modulo rescaling, X ~ aX. The six-dimensional coordinates X 
modulo rescaling can be interpreted, in a standard way, as parametrizing anti-de-Sitter space. 
The null condition X'^ = then gives the usual realization of (conformally compactified) four- 
dimensional Minkowski space with the boundary of AdS space (see, for instance, [35] for a 
pedagogical presentation). The conformal field theory and its fields are functions on this 
four-dimensional boundary. Thus we are simply using convenient coordinates to describe 
Minkowski spacetime. 

A procedure for describing general spin fields in this language was given in [34]. We will 
be using a slight variation of this method. The simplest example is a scalar field with mass 
dimension 1, which becomes a function (p{X) with homogeneity degree —linX. Its two-point 
function, up to normalization, is then the only invariant with the correct weight 

iHXmX)) = (2.2) 

where {XY) is the six-dimensional dot product. As in [34], we will need to lift 2-component 
4-dimensional Weyl spinors ipa, to 4-component 6-dimensional Weyl spinors ipa- We will do 
this in a slightly different way, though. We choose to remove the extra two components 
through the equivalence relation ip ip + X(, where the second term involves gamma matrix 
multiplication. Because X'^ = 0, the space of spinors of the form X(^ has dimensionality two 
as required. In general a field with s spinor indices and dimensionality d will be mapped 
to a field with homogeneity degree —{d+^s). Thus, a canonical fermion field will have the 
propagator^ 

(^.(x)^^(y)) = (2.3) 

^Our fermions are related simply to those of Weinberg [34], which satisfy the constraint X-^woinborg = 0: 

^Weinberg ~ XtpuB- 
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Our notation is closely related to the "momentum twistors" of Hodges [36] . If X is a null 
six-vector, the antisymmetric matrix X"'' has rank two and can be written 

for some four-component "momentum twistors" X" , X"" . These are the two solutions to 
the Weyl equation Xip"- = 0. 

This is the usual twistor correspondence between a point in spacetime and the line 
in momentum twistor space. (This is the standard twistor correspondence. We 
use the name "momentum twistor" , as opposed to simply twistors, only as a reminder that the 
spacetime we are working with, on which the Wilson loop is defined, corresponds physically 
to momentum space in the scattering amplitudes side.) 

The momentum twistor variables are the most efficient way to parametrize a null polygon 
contour, or, equivalently, the data of a massless planar scattering S-matrix element. Let us 
quickly review the mapping, due to Hodges [36]. The external momenta pi are related to the 
dual coordinates Xi through 

pr=xr-xf_,, (2.5) 

where Xn+i is periodically identified with xi. Each point Xi is a-priori associated with 
two momentum twistors, but the null condition (xj — Xj-i)^ = means that Xi and Xi-i 
share a momentum twistor; call it Zj. Therefore the null polygon is described by n free 
momentum twistors Zj, and its corners are located at Xi = Zi A Zi^\. The same can be done 
supersymmetrically. Writing p"" = A"Af for spinors z, and using Grassmann variables to 
describe the external polarizations, one introduces 0's through 

A?^.^ = ^r^-e1. (2.6) 

The momentum supertwistors are then, explicitly, (Zi^Xi) = (Aj, (AjXj), {Xi9i)). 

The dot product between two 6-vectors is precisely the fully antisymmetric contraction 
of the four momentum twistors (XY) = {X^^ xj^^Ya^^Ya"^^), which explains the bracket 
notation. We will repeatedly write strings of six-vectors acting on a momentum twistor, for 
instance objets like {0XYWe'). UY = F^^) A T^^), this is to be interpeted as 

{exYwe') = (0xy(i))(y(2)i^0') _ {exY'^^^){Y'^^'>we'). (2.7) 

These products obey the Clifford identity {{X,Y}) = —{XY). The remaining rules needed 
to convert a general space-time expression to six-dimensional notation are: {x — y)'^ — )• 
{XY) / {X I) {Y I) , and ^ J d'^x{XI)^ — )■ f-^, where I is the infinity point. The spacetime 
integral is normalized so that 



Most integrals ever needed can be derived from this basic one or its derivatives, thanks to 
Feynman parameterization — a concrete example will be given in subsection 4.3. 
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Derivatives d/dX are well-defined modulo amounts proportional to X. For instance, the 
following is the standard fermion kinetic term in the action 



^2 , V'^^aaV'" ^ / "P^Q^^ + I ^XAXi^. (2.9) 



It is verified that this is invariant under the substitution ifj — t- C,X. Integration by parts, using 
the general identity fxi-^wc ~ ^)-^ ~ ^' verified to relate this expression to a similar 

one with ^lJ and -0 exchanged, showing that the unphysical ^ polarizations also decouple. One 
has to use the homogeneity condition {X-^)il; = —2ip. 

2.2 Chalmers- Siegel action and gauge-fixing 

The Chalmers-Siegel formulation of gauge theory reorganizes perturbation theory as an ex- 
pansion around self-dual Yang-Mills. This arises naturally in the context of supersymmetric 
Wilson loops, where supersymmetric Wilson loops computed within the self-dual theory are 
found to be dual to tree-level N^MHV tree-level scattering amplitudes [29, 30]. This will 
also prove quite useful for the two-loop computation we present below. The reorganization 
proceeds by rewriting the action as 

^9ym J ^9ym J 



47r2 



j d'xF^pG^^ - j d'xG^pG^P. (2.10) 

The first line corresponds to the addition of a (imaginary) ^-term, which has no effect in 
perturbation theory, while the second line introduces an auxiliary self-dual field G. It is 
natural to perform similar manipulations with the gauge-fixing term in Feynman gauge: 

These two equations are unified if we upgrade the spin-1 fields to general two-by-two matrices: 



Go.p ^ + le^^A, F^p ^ da^^Ap^ + [A^^, Ap% (2.12) 



The propagator between A and the new G is simply 



{AUx)Gp,{y)) = ^ (2-13) 

This propagator is not conformally invariant. Following the rules of the six-dimensional 
notation, G is mapped to a weight —3 bispinor and the propagator becomes 

{A,dx^{X){9Ge'){Y)) = f^y^f^YI) ^^-^^^ 
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where we have contracted G with two fiducial momentum twistors 6 and 9' lying on the line 
Y . The action of conformal transformations on the gauge-fixing generates a four-dimensional 
family of gauges, labeled by different points I's. (A five-dimensional extension where 7^ 
might also be consistent, although we will not consider it.) This freedom to choose I will be 
exploited below. 

2.3 Review of loop integrands 

Let us conclude these considerations with a simple example illustrating the formalism at 
work: the 1-loop correction to a bosonic Wilson loop (which is dual to a MHV amplitude). 
To get the order g^-^ correction to the Wilson loop, we need to bring down one power of the 
interaction Lagrangian and take its correlation function with the Wilson loop. First we need 
the correlation function between the Wilson loop and a single field strength at the point Z.^ 
This receives one term per edge. Parametrizing edge i using X = Xi-i + TxXi, this gives 

The integral gives simply 

^ (2.16) 



{Zi-U){Zii+l)' 

while the prefactor is simply {i—lii+19){iI9'). Therefore that correlation function is 

The sum can be verified to be independent of / using the four-term identity 

{i-Ui+19){Ii9') {i-lii+l9){I'i9') {9n+lII'9') {9i-liir9') 



{Zi-li){Zii+l){ZI) {Zi-li){Zn+l){Z'I) {Zii+l){ZI){Zr) {Zi-li){ZI){ZI') 

(2.18) 

which transforms it to a telescopic sum. To obtain the loop integrand we substitute this into 
the interaction Lagrangian (the order gy^i part of the Chalmers-Siegel Lagrangian), which 
reads 

Sint = -9\ GabZ-'Z^'^Gcd + .... (2.19) 

where = ^^^^ and the dots denote Yukawa and scalar self-interaction terms. After using 
the simple identity {Zij){Z{Ii) n (Ij)) = {Z{i) Pi {j)){Iij){ZI) one obtains, directly^ 

,i-ioop^_ V f {Z{i-ln+l)nij-ljj+l)){Iij) 

,J^^^Jz{Zi-li){Zn+l){Zj-lj){Zjj+l){Ziy ^ 

This is exactly the 1-loop MHV integrand in the form found in [23]. The integrals will be 
discussed below. The parity even part of this expression reproduces the more familiar sum 
over "two-mass-easy" box integrals [37] (this is verified in detail in [31]). 

^In this paper we will sometimes use the letter Z without a subscript denotes a spacetime point; we hope 
this will cause no confusion with the momentum twistors Zi. 

*The intersection symbol is defined as {ij{abc) PI (def)) = {iahc){jdef) — {jabd) (idef) 
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2.4 Symbols 

Transcendental functions can typically be written in many equivalent ways, which usually 
makes it difficult to decide whether two expressions are equal or not. A spectacular recent 
example is the equality between the compact formula given in [1] and the original 17-page 
analytic formula obtained in [2], for the hexagon 2-loop remainder function. 

The "symbol" of a transcendental function is a piece of invariant data which trivializes 
all identities between transcendental functions. The symbol is uniquely defined and if two 
expressions have different symbols, one can immediately conclude that they differ. The symbol 
is lossy, however, and from the equality of two symbols, one can generally only infer that they 
agree up to lower functional transcendentality, and choice of Riemann sheet. 

This subsection is meant as a lightning introduction to symbols and as a refreshment for 
the reader's memory. Nice presentations can be found in [1, 15]. The transcendental functions 
of interest to us are (generalized) polylogarithms and are iterated integrals of the form 

F{x)= I d\ogX{xi)d\ogY{x2)d\ozZ{x^). (2.21) 

J Xq<X\<X2<Xs<X 

This particular integral would correspond to a transcendentality-3 function of x, e.g. some 
combination of Lis, Li2 log or log^ with various x-dependent arguments. The variable x and 
the integration path in general live in some multidimensional manifold; for us it is useful to 
visualize a point x as a configuration of momentum twistors, or, equivalently, a collection 
of cross-ratios. The symbol is simply the integrand of such an iterated integral, written as 
log X (S> log Y (S> log Z, and abbreviated to 

S[F{x)]=X0Y0Z. (2.22) 

The integrals of interest are such that they are contour-independent, expect for monodromy 
when the contour crosses a pole. The upshot is that all identities between transcendental 
function descend to trivial identities at the level of the integrand, and are mapped to simple 
algebraic identities on the symbol. The only manipulations allowed are multilinearity, e.g. 

X (g)YW (g) Z = X 0Y Z + X (E)W Z, (2.23) 

together with ^[Constant] = 0. Frequent examples include: 5[log x log y] = x®y + x and 
5[Li2(l-x)] = -x®{l-x). 

3 Supersymmetric Wilson loops 
3.1 Chiral Wilson loops revisited 

It is well-known that supersymmetric gauge theories can be interpreted as gauge theories 
in superspace. Unbeknownst to the author of [30] at the time of that work, is the equally 
well-known fact that N = A Yang-Mills can be interpreted as a gauge theory in a superspace 
(x, d^, 9\), where A = 1 ... 4 are SU(4)ij indices [38, 39]. In that case, the superconnection 
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is defined only on-shell. In other words, given any on-shell solutions to the field equation, 
one can construct a superconnection which is supersymmetry covariant. Actually, it descends 
from a superconnection in ten-dimensional M = 1 super Yang-Mills. 

A chiral supersymmetric Wilson loop, for null polygon contours, was constructed in [30] 
as follows. The bosonic Wilson line along edge i is translated into superspace: 

where c = ( 2^^^ Y^'^- This result is not quite yet supersymmetry covariant, e.g. not invariant 
under cq'\ — Ylii Ka-J^- The reason is that in attempting to verify invariance one needs to 
commute q across the exponential, and one picks up an anticommutator 

{g^,g|}A^ = e„/3Z)^</.^^. (3.2) 

In [30] this non-invariance by a total derivative was canceled by inserting operators at the 
cusps, and the choice A=k = Aj_i was made.^ 

Since the defining property of the object we are interested in is its supersymmetry covari- 
ance, and since an object with this property is obviously unique, the result is guaranteed to 
be equivalent to the standard = 4 superconnection with 6 set to zero in it, up to equations 
of motion and gauge transformations. 

We can even construct explicitly the super-gauge transformation which relates the two 
(a similar computation was carried out in [40]). From the integrability conditions obeyed 
by the standard superconnection, F^^^^-^ = and P^^p) ~ ^' ^^^^^ (o^/3) means symmetrical 
part, one concludes [29] that the ^-path on top of edge i can be deformed by any amount 
proportional to Aj. If we write the 0-path as 

e{t) = (1 - t)6i_i + tOi 

-t\i_i{u+l) + {1 - t)\i+i{i-li) f,,.Xi-ih , .Xi+ih , , /,NXjAi(i-li+l) 
= Xi /. T\/--ii\ + / ^ /■ T\ - gw /■• , i\ + W) /■ T\/- - n\ 

where f{t) = (1 — t), g{t) = t and h{t) = are the canonical choices, what we are saying 
is that different functions /(t), g{t) and h{t) are simply different gauge choices. The gauge 
taken in [30] corresponds to step-functions f{t) = 9{e — t) and g{t) = 6{t — 1 -|- e) with 
e very small. This made manifest the fact that Xi-i Xi+i have no coupling inside to 
edge i. (In addition, the choice h{t) = (1 — t)9{t — e) was made in [30], which removed the 
Aj+i component. In retrospect, the choice h{t) = would have been much better, as it is 
conformally invariant.) 



""From this viewpoint we can propose a better starting point for a nonperturbative definition: conjugate 
with 

e <"-i"+i> (3.3) 

where Q is the 4-component superconformal generator. Q will be slightly broken at the cusps by infrared 
regularization, but this can usually be controlled. At least naively this definition removes cusp terms. The 
author thanks J. Maldacena for related discussions. 
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We conclude this discussion by reporting the beautifuhy simple expression we obtain for 
the abelian superconnection in 6- dimensional notation, 



A = {Adz) + ii^edz) + {9de—){(j) + 9i; + ece). 



(3.4) 



up to signs and numerical factors {O"" is a 4-spinor with Z9 = 0). 
3.2 A conjecture 

It is now easy to state our proposal, or conjecture: 

Conjecture There exists a generalization of planar scattering amplitudes, An{xi,6i,9i), 
such that 



plus similar equations for the dual special superconformal generators, for any infrared-safe 
quantity, such as the remainder function or the ratio of two A^s with different MHV degree 
but at the same kinematic point. Furthermore, this object should be equal to the expectation 
value of the supersymmetric Wilson loop of [39] on a null-polygonal Wilson loop contour. 

To be fully precise, we should state what the "null" condition is, because it enters the data 
in An- The condition should leave 3 bosonic degrees of freedom per edge as usual, plus four 
chiral fermions and four antichiral fermions, and there can be only one such superconformally 
invariant condition. Tentatively, we would write it as pf = 0, Piaa{Gi — = 0, Piaai^i — 

9i-i)°' = where pi = (xj — -|- 9i9i^i — 9i^i9i). However we do not want to embark 
on checking these conditions here, because they will not be needed in this work, and present 
them only as a tentative guess. 

A more important question, in our view, is to try to write these Ward identities in terms 
of unconstrained variables such as momentum twistors. To first order in the new variables 
the following algebra is easily deduced: 




(3.5b) 



(3.5a) 




(3.6b) 



(3.6a) 



where a = 1 ... 4. This is the algebra we will be working with in this paper. This does 
not look very parity symmetric, but is simply an artifact of working with chiral variables. 



Unfortunately, we do not know yet a nice extension of this equation to higher orders in 

In terms of these variables, the first few terms of the super- Wilson loop are, in the 
normalization of [30], 

A = dxaaiA'^" + + sV^" + •••), (3-7) 
where 6i = {xi^i+i — or, in 6d notation with X = Xi^i + TxXi, 

A = {dXA) + Xi{i-^n+'i^i^)dTx + xS'4')d'Tx + • • • (3.8) 

The dots also contain a term F06 which would contribute to the same order in XX- Fortu- 
nately, we will manage to avoid needing this term. 

The reader will notice that we have also not specified what the quantum definition of the 
Wilson loop should be (e.g., including its proper regularization.^) Only its bare expression is 
known. Obviously, the defining property of the object we are interested is the Ward identities 
it should obey, and part of the conjecture is that a regularization preserving this property 
exists. 

We mention at least one, in all appearances natural, quantum interpretation of the clas- 
sical Wilson loop which leads to e/e-type violations of supersymmetry [40]. In other words, 
regularization appears to be not entirely straightforward. A regularization of the (square of) 
chiral Wilson loops based on correlation functions of BPS operators which corrects for these 
problems was proposed thereafter in [32, 33]. We will sidestep these questions in this paper 
by asking simply whether a consistent answer obeying the Ward identities exists. Normally, 
this is the right question to ask whenever one is worried about anomalies — we know of no 
example in quantum field theory where a consistent answer exists and yet there is an anomaly. 

3.3 MHV case 

We now concentrate on the MHV remainder function. With no loss of generality we can 
expand 

rT"" = RT^\^,^=o + E liJ^u+l-f ''^ + (adjacent terms) + 0{{xx?). (3.9) 

non-adjacent 

The "non-adjacent" terms cover all i ^ j, For adjacent i and j, we define Cij = 0. 

Let us work out the consequences of the Ward identities. From {i—lii+lQ) Rn = 
(no summation over i), we deduce the sum rule 



Y^Cij = 0. (3.10) 



''One promising possibility, suggested to the author by Beisert, Huang, Vergu and Skinner, is to introduce 
two dual supertwistors Z and W subject to the superconformal constraints Zi ■ Wi =0, Zi ■ W,;_i — and 
Zi ■ Wi+i — 0. The constraints can be used to solve for the bosonic W's. The symmetry generators are then 
simply Zd/dZ + Wd/dW. 

This does not mean an off-shell definition, which obviously would be too much to ask for in A/" = 4, 
but which is not required either. It should be perfectly legal to use equations of motion inside a correlation 
function provided 5^{x) terms can be neglected, as is generally the case for Wilson loops. 
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Similarly, there is the parity conjugate sum rule 



^Cij=0. (3.11) 

i 

From * Q)Rn = 0, we deduce that 

"'^'^=° - z^^*'^(.-Hi+ij)(i-iu+i*) ^^-^'^ 



which shows that all terms linear in x ^md x ^'^^ determined by the Cij^s. For instance, a 
similar equation using * Q) would determine the coefficient of XiXi- The right-hand 

side is independent of *, thanks to the sum rule (3.10). 
From -^QRn = 0, we obtain that 

3 i ^ J 

This shows that all bosonic first derivatives of the remainder function are determined by the 
order xx terms. Eq. (3.12) allows us to express the right-hand side in terms of the Cjj's solely. 
Let us contract both sides of the previous equation with dZ^ and sum over j. Collecting the 
terms involving Cj.j gives 

= Cij(ilog(i-lM+lj). (3.14) 
We have thus obtained a simple result for the total differential: 

dRn\x,x=o = ^ Cijdlog{i-ln+lj). (3.15) 

This is the main result of the present subsection. The left-hand side is the total differential of 
the bosonic remainder function we are interested in. The coefficients Cj.j are defined by the 
expansion (3.9). As explained in the previous subsection, they may be computed by turning 
on fermions on specific edges of the Wilson loop. 

Let us try to compute the Cjj's, say C2,i, at one-loop. We have a ^p propagating along 
edge 2 and a ip propagating along edge i (see Figure 1), joined by a propagator: 



(2i)(i2) 
{XYy 



C2,i= I drxdry^^^—y = log n2,i--i,i,i (3.16) 



where X = {Z\ — txZ^) A Z2 and Y = {Zi^i — TyZi^i) A Zf, we define cross-ratios as 

_ {xa,){x,xi) 
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i-] 



Figure 1. The diagram giving the coefficient function C2.i at f-loop. 

We note three facts. First, the sum rule Eq. (3.10) is almost obeyed, the sum over i being 
telescopic. However there are some uncanceled boundary terms. Second, C2,4 and C2,n are 
ill-defined, plagued with infrared divergences. Third, leaving aside these boundary issues, 
the general term in Eq. (3.15) is easily verified to perfectly match with the derivative of the 
1-loop MHV amplitude [37]. 

The slight mismatch here with boundary terms was to be fully expected: the naive Ward 
identities should be modified to account for the infrared divergences, probably in a very similar 
way to the bosonic case [41] . We note that there exists a canonical infrared regulator of the 
theory, the Coulomb branch, in which these anomalies are removed by an appropriate action 
on the moduli space parameters [42]. It would be tantalizing to interpret the x variables as 
some sort of supersymmetrization of that regulator. For the moment, we concentrate on the 
remainder function, which is infrared finite and for which the naive Ward identities are fully 
justified. 

3.4 Leading singularities 

The fact that MHV amplitudes turn out to be sums of pure transcendental functions (meaning, 
transcendental functions with rational prefactor equal to 1), is generally understood from the 
fact that the only Yangian invariant with no fermion is a constant. The reasoning takes its 
roots in generalized unitarity, and is based on the idea (highly plausible, but for which there 
exists no general proof at the moment) that the coefficients of the transcendental functions 
should correspond to "leading singularities" of loop integrands, that is, by the loop integrals 
performed on closed compact cycles (see, for instance, [43]). These leading singularities have 
been proved to be Yangian invariant [23, 44]. 

Therefore, a key question to ask about the Cjj's is, what are their leading singularities? 
That is, what can arise if one performs the integrals contributing to a Cj j on a closed compact 
cycle? 

In the case of x = 0, all invariants are known to be closed contour integrals in the 
Grassmanian [45-47]. In the MHV case, this meant the constant '1'. The computation in the 
next section strongly supports the idea that the Cjj 's enjoy similar properties - their leading 
singularities at two- loops are the same as at one- loop, and are all '1'! (The fact that one 
finds c-numbers, as opposed to rational functions, is the main statement. However, we do not 
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even find any g-) This immediately implies that the last entry of the symbol, at both one- 
and two- loops, is of the form {i—lii+lj). 

It would be very interesting to see if this property is general - whether leading singularities 
in the presence of x's at N'^MHV for a given number of external legs form a finite set to all 
loop orders -, and whether in the MHV cases these leading singularities are indeed saturated 
at 1-loop.*. 

4 Two-loop MHV remainder function 

We wish to use Eq. (3.15) to determine the derivative of the two-loop n-point remainder 
function. 

4.1 Outline of method 

There exists of course many different and equally valid ways to calculate the Feynman dia- 
grams which arise. We have employed a certain systematic method, which we now describe. 
The advantage of the method is that it automatically leads to relatively nice one-dimensional 
integrals, with manifest transcendentality degree 3. In our view, this by far compensates for 
the somewhat more lengthy algebraic manipulations the method entails. We also believe it 
likely that the method has a generalization to higher loop orders. 
The steps of the method are the following: 

• We fix the parameters tx and ry that describe the insertion points of the fermions, and 
compute the integrand associated with a loop insertion point Z using the method of 
Lagrangian insertion. This gives C2,i as a rational integral over 6 variables: tx, ty and 
the space-time point Z. 

• We perform the integration over ry. The motivation for doing ry first is the following. 
The fact that the Lagrangian is chiral means that it has different types of singularities 
near X and Y. For instance, the OPE of the chiral Lagrangian with ip and ip dictates 
double poles in (ZY)'^ but only single poles in (ZX). Pleasingly, the integration over 
ry will simply convert the double poles to single poles, leaving a rational expression. 
At this stage we thus have a five-dimensional rational integral with single poles. 

• We perform the four integrations over Z. This has the structure of a standard loop 
integral and produces dilogarithms. In fact, all the integrals will be known pentagons 
and scalar boxes. 

We see that the problem is reduced to a one-dimensional integral over dilogarithms, 
without the need for further "art" or trickery. 

We found the following technical tricks to be very useful: 

*In the first version of tliis contribution, a symmetry argument was suggested based on an (incorrect) 
classification of "invariants of iiomogeneous degree in x and x" ■ Invariants witli this property do not actually 
not exist when the full form of the superconformal algebra is used. For this reason this discussion has been 
removed in second version. 
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We split the computation into four individually gauge-invariant and conformally invari- 
ant subsets. This involves adding and subtracting certain "geodesies in AdS" edges to 
the polygon. While this may sound a bit unfamiliar, the idea in picture is very simple: 





C 

On the left are the graphs contributing to the logarithm of the amplitude: The gluons 
couple to all edges of the two shaded regions and the Lagrangian is inserted somewhere 
along the wigggly line. (The figure includes a graph in which the two gluons are attached 
to the fermion propagator.) On the right, each polygon is split into two, one part outside 
the box (l,2,i— and one part inside, leading to four terms. The reason we can do 
this so simply is that the gauge sector is essentially abelian. 

Making sense of these figures requires specifying a path connecting X2 to Xi-i. Since 
these points are not null-separated, no path in spacetime will do this in a conformal 
invariant fashion. We want to preserve conformal invariance. The way out is to use a 
path that goes through a fifth-dimensional AdS space. Let us explain what this means. 
It is here that the six-dimensional notation introduced above pays off. The gauge link 
from two null-separated points Xi and X2 was written in section 2 as 



{6XiX2ie') 

{ZXi){ZX2){ZI) 



(4.1) 



When we verified independence with respect to / there, nowhere did we need that Xi 
and X2 where null separated (nor even that they were points inside spacetime, Xf = 0). 
Thus for general Xi and X2, we simply define the gauge- link along an AdS geodesic to 
be that expression. This provides a gauge-invariant closure of the polygon. 

We render the top and bottom polygons finite in the simplest possible way, by moving 
their endpoints slightly away from the corners of the box: 



In equations, 

Z2AZ3 Z2l^Z^ = (Z2-eiZ4) AZ3, Zi_iAZj Zj_iAZj/ = A (Zj -e2^^-2)• 
(4.2) 

The lower polygon is to be regulated in the same way. C^^^ is trivially finite as the 
regulators are taken away, but C*^^^ and C^'^'^ diverge logarithmically with ei and €2- 
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We will see, however, that these divergence can be subtracted and renormalized in a 
canonical and conformal invariant fashion. Thus all infrared divergences will be shoved 
into CW. 

• We would still have to regulate the graphs contributing to C^'^\ should we want to 
calculate them. And we would have to carefully carry out the subtraction of the two- 
loop correction to the cusp anomalous dimension that appear in the BDS Ansatz [5], 
following the definition of the remainder function. However, by construction, C^'^j will 
be a conformal invariant function of a single cross-ratio Uij-ij^i-i. So, instead of 
computing it directly, we can use supersymmetry to obtain it from what we already 
have, through the sum rule (3.10). This way, the whole computation is phrased directly 
in terms of the remainder function. 

4.2 Contribution C^^^ 

To obtain the integrand we begin by writing down the field strength induced by the top and 
bottom polygons 

(nrf)') - {j-ijj+ie){jm _ {e23i-iii9') 

(f^G^y Jtop - iZj-lj){Zjj + l){ZI) {Z23){Z^-l^){ZIy ^^'"^ 

_ v^' {k-ikk+w){ki0') {euii+iie') 

jbottom - ^Z_^^ {Zk-lk){Zkk+l){ZI) ^ {Z12){Zn+l){ZI) ■ ^ ' 

The spacetime point Z will be the loop integration point and 6, 9' are two arbitrary 4-spinors 
associated with that point. We see two sorts of terms. The terms inside the summation 
come from the null edges and are as in Eq. (2.17). The extra term comes from the "AdS 
edge" mentioned above, and ensures that the field strength of each polygon is gauge invariant 
(independent of /). We should remark that there is some freedom regarding the ordering of 
the first two factors in these terms, e.g. XdX versus —dXX, which did not arise in the case 
of a null edge. Any choice would be just as much gauge-invariant. The choice we have taken 
turns out to ensure the finiteness of the integrals, as will be fully apparent shortly. It would 
be nice to understand better what is the special physical meaning of this choice. 

The simplest way to compute the product GtopGbottom is to make the gauge choice 
/ = (ii+l) in both polygons. This kills two edges in the bottom polygon and yields the 
integrand 

V V ha- 7-1 r k-l k) - V / {Zi23^)nim^-ln+lk) 

^i^Jz{Z23){Z^-U){Zn+l){Zk-lk){Zkk+ly ^^''^ 



This multiply the fermion exchange diagram computed in Eq. (3.16). The integrand is easily 
verified to be finite, and can thus be evaluated without regularization. Term by term it 
contains collinear singularities on edges i and i+1, however. What appears to be the simplest 
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way to proceed is to regulate the individual terms by simply moving (ii+l) — )• (ii+l) + el' . 
The basic integrals that we need are 



h{i,j;i) 



{Z{i-iii+i)nij-ijj+i)){iij) 

z {Zi-li) {Zii+l){Zk-lk) {Zkk+l) {ZI) 

Li2(l - U/,l,2,j) - Li2(l - U^i^2,i-l) + Li2(l - W2,j_l,j,l) 

+Li2(l - - Li2(l - nx,i,i-i,2) + logiix,i,2,i logux,i,i-i,i (4.6) 



and 



r2me' 



2Li2(l -t;) +logwlogw + C'(e), (4.7) 



where v = Vij-ij^i^i and w = Ui'^i^ijj^i' oc e^. Indeed, using the Schouten identity 

(Z(23i) n (fc))(2A;) = (Z(i) n (^))(23«A;) + {{Zk-lk){kk+l(i) n (i23)) - (A;-l ^ k+1)), (4.8) 
everything is reduced to these two basic integrals. The final result is written compactly 



i-l n+1 

Ci]} =logU2,i-l,i,l xY^Yl 
j=2 k=i 



Li2(l - n^- fc,(fc_i),(j+i)) + log log 2 ^'^ 



(4.9) 



if we use the notation that = X2 when j = and X(fc_i) = xi when k = i. This is 

quite similar to the pairing of two Wilson loops in [15], except now we have two edges in AdS. 
(There is a slight abuse of notation here - some of the dilogarithms and logarithms in the sum 
are divergent. They can be removed using the identity Li2(l — u) = — Li2(l — 1/u) — ^ log^ u, 
at the cost of making the expression look less symmetrical.) 

4.3 Contributions C^^) and C(3) 

To find this integrand, we now need the three point function for a field strength and two 
fermions. The Feynman rules from Eq. (2.9) gives this as 

_ f {eu2){u(i)niie')) 

'^'^ Ju {UX)^{UY)^{UZ)^{ZI)' ^ ^ ^ 

The integration over the spacetime point U can be done using the standard Feynman pa- 
rameterization procedure. The six-dimensional notation offfers a particularly efficient way to 
achieve this; let us go through the various step in detail (see also [48-50]). One first combines 
denominators using Feynman parameters 

r(„) f „„^^ '"r^;--!. (411) 



(f/Xi)((7X2)---{(7X„) ' 'J {U{o.iXi+a2X2 + ...+o.„X„)y 

The measure [d^^^ai ■ ■ ■ an] is the projective measure, defined explicitly by setting one of the 
coordinates to 1 and integrating over the n— 1 others. It does not matter which variable is 
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set to 1. In the present case, the loop integral that we need is the second derivative of the 
basic integral in Eq. (2.8): 

f {CDY){CDZ) _ {XY){XZ) _ {YZ) 

Thus Eq. (4.10) becomes 

[d?abc]abc f6{e{bY)2){{aX + cZ) (J) n {W')) 2(^2 (J) n {16')) 



Vx. 



{ZI) \ {ab{XY) + ac{XZ) + bc{Y Z))"^ {ab{XY) + ac{XZ) + bc(Y Z)Y 

1 ( {QY2){x{i)r\{w')) {eY2){z{i)n{ie')) \ 

{XY){XZ){YZ){ZI) [ {XY) + (YZ) ^ ) 

{eY2){i9') ^ miexYie') ^^^^^^ 



{XY){XZ){YZ)^ {XY)^{XZ){YZ){ZI)' 

In the last line we have used Schouten identities to rearrange the result to a nicer form. 

As anticipated, we find a double pole in (YZ) but only single poles in (XZ). We are 
going to perform the ry integral first. Before doing that, we should combine this term with 
the segment from Y to Xi^i, since this side of the square also has a nontrivial dependence 
on Y. Parametrizing Y = (Zj_i — ryZj+i) A Zi, this gives 

' / {i2){i-ln+19){zI9' 
dry vxy,z + 



{ZY){Zii+l){ZI){XY)'^ ^ 
{ie){i-lie'2) {z2){en+lXie') 



(Xi-li) (XZ) (Zi-li) (Zii+l) (Xi-li) (Xii+l) (XZ) (Zii+l) (ZI) 



(4.14) 



The fact that this is rational lies at the heart of our approach. The remaining two sides of 
the square are then obtained as in the previous section. After adding them, the / dependence 
disappears in a simple way, yielding finally 

inrf)') = {i9){i-ue'2) {i2){912ii+lX9') 

{ Jbottom square (^Xi-li) {X Z) {Zi-li) {Zii+l) {Xii-1) {Xii+l) {X Z) {Z12) {Zii+l) ' 

(4.15) 

We are nearly done with the computation of the integrand. The final step is to contract 
this result with Gtop which was obtained in the previous subsection. Following the adopted 
strategy, this is to be done with (23) — )• (23)' and {i—li) — ?• {i—li)' in order to regulate 
collinear divergences. Because everything is then safely finite, we can simplify our life by 
making the gauge choice / = {i—li)' inside Gtopj which removes a couple of terms. A short 
computation using Schouten identities then reduces everything to the known integrals listed 
previously: 



(4.16) 
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where 



J = I dxlog^^^(/5(^;2',3;i-l,i)-/2-^'(X,2';i-l',i)) 

+ I dxlog l^^'^^^'l [h{i; 1,X; 2^ 3) + h{i; 1, X; i-2, i-l') + l'"''' {X, 2'; i) 

-/2^^'(l,2',i-l',i)V (4.17) 



This completes the computation of the integrand associated with the two-loop derivative C^^^. 
The next step is to renormalize this expression. When the corners 2' and i' of the polygon 
approach their physical values 2 and i, J diverges like 

/ ^ /"\ (Xi-li) , 
\og{i-li ii+l) J dxlog -^^^^ log ni,2,x,i-i 

+ log(122'3) r dxlog iy'r^'l logn,_i.i,x.3. (4.18) 

These are all the divergences. A useful fact, easily seen by integration by parts, is that 
these two integrals are equal, e.g., the logarithms combine into log((i— li'ii-|-l)(122'3)). 
This was guaranteed by symmetry. Importantly, this means that the divergence can be 
renormalized in a simple and conformal invariant fashion, by replacing this logarithm with 
log((12z—li) (23^+1)), which has the same little group weights. We have thus succeeded in 
making this subset of the computation conformally invariant on its own. 

We can do better. As discussed in the preamble, in the case i = 5 where the top polygon 
is a triangle, conformal symmetry ensures that C2 5 depends only on ^1,4,5,2- In general, it 
makes sense to improve our renormalization prescription by adding a counter-term depending 
only on U2,i-i,j,i, which means moving things between C^^-* and C^'^\ This allows us to choose 
a counter-term such that the triangle vanishes. Although we will not need it, we record for 
completeness the expression for the counter-term 

/■ , , (Xi-li) , 

Jet. = - J dxlog -^^^^ log Ui,2,X,i-l log Ui,2',j-l',i 

+ J dxlog (2Li2(l - iij,i,2,i-i) - 2Li2(l - Mi,x,2,j-i) - l^h{1 - ^ij,i,x,j-i) 

-hLi2(l -Ui^i^x,3) -Li2(l -ni_i,i,x,3) -h log ni_i_i,x,3 logtii,2,x,i-i) 
+ J dxlog -^^^^ (Li2(l - nx,i,i-i,2) + logtix,j,j-i,2 log ux,3,i-i,i) • (4.19) 

Apart from the first line this depends only on the cross ratio U2,j-i,j,i-^ What is more impor- 

^ This does not depend on the point 3 because it enters only through the combination {X34)/{1234). In 
our parametrization of X , this is equal to 1. 
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tant, and can be obtained directly without working out Jcx.^ is the renormalized expression: 

(Xi-li) 



J + J, 



c.t. 



/ 



(Li2(l - nj^i^x.s) - Li2(l - Mi-i,i,x,3) - Li2(l - 



+Li2(l - ni_i,i,x,i-2) - log'Ui,3,i_2,i logUi,i,x, 



+ / dxlog ^^3^^^^ (Li2(l - ^t2,i,i-l,3) - Li2(l - 1iX,i,j-l,3)) • 



(4.20) 



This vanishes for i = 5. Substituting this for J in Eq. (4.16) and rearranging a telescopic 
sum, we conclude that 



C. 



(2) 
2,j 



*-2 .I 



^ / dxlo^ 



,=4-1 
«-2 /.2 



+ V / dxlo; 



{Xi-\i) 
{Xi-li) 



h{X;i-l,i;j-l,j) 



(Li2(l - nj_i,i,xj) + Li2(l - Uijj-i^i) + Li2(l - 



{Xii+1 

-Li2(l - Ui_i,i,xj-i) - Li2(l - Uijj^i^x) + logiii_i,i,x,i log 



{XM) 

4.4 Contribution C^^^ 



(Li2(l - U2,i,i-l,3) - Li2(l - Ux,i,i-1,3)) ■ 



(4.21) 



As explained above, we don't actually need to compute anything new to obtain C^^^. By 
construction, C^'^j' depends on a single cross-ratio and can be determined through the sum 
rule (3.10). Let us denote this cross-ratio as n = Uij-ij^i-i. From considering the sum 
rule for the pentagon, where the remainder function and the contributions C^^\ . . . , C^^-* all 
vanish, we conclude that C^'^\0) = 0. In the hexagon case, the nonzero terms in the sum rule 
give 

Ci%) = -Ci]l - Ci$ - Ci% (4.22) 

It is a nontrivial consistency check that the right-hand side depends only on the required 
cross-ratio; if that were not the case, the proposal would be wrong. 
In the hexagon case, Eq. (4.9) evaluates to 



C^^l = logU X (Li2(l -U3,5,6,2) + Li2(l -U4,6,l,3) + log ^3,5,6,2 log U4,6,l,3 " ^ 



(4.23) 



Representing log u as dx log j^^gy , this can be added directly to the integral representation 
(4.21) of and C^^l We find: 



dxlog l^ll^ (Li2(l - Ux,5A,2) + Li2(l - UXfi,5,2) - Li2(l - liX,6,4,2)) 

+dxlog )|||( (log Ux,5,6,2 log nx,3,4,l) 

1X61) 
(X56> 



+dxlog (Li2(l - Ux,4,5,l) + Li2(l - ■wx,3,4,i) - Li2(l - Ux,3,5,l)) 



.(4.24) 
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The three hnes are actually all equal in magnitude, but alternating in sign — this can be 
proved by integrating by part in the first or third line. So the total is minus the second line. 
It can be expressed in terms of the required cross-ratio by a simple rescaling of tx '■ 

Ci%) = f^~^ d(\og ^) log(l + r) log(l + ^) 

= _2Li3(l--)-Li2(l--)logu-^log3tx+^logn. (4.25) 
u nob 

The same expression (which is valid for u > 0), with the appropriate cross-ratio, gives all 
4.5 Method for computing the symbol 

We still have to integrate Eq. (4.21). The first step is to compute its symbol. It is actually 
possible to do so without performing any integration. We are not aware of any standard 
algorithm for doing so, but we will describe the method we have employed. The method is 
based on computing discontinuities across branch cuts, and comparing with the discontinuities 
of the iterated integral (2.21). The leading transcendentality branch cuts of the latter end 
at the zeros and poles of the leftmost entry of the symbol. These discontuinuities themselves 
have branch points, which are at the zeros and poles of the second entry of the symbol. Thus, 
by computing discontinuities of discontinuities and so on, one can read off the symbol. 

We hope to elaborate elsewhere about the algorithm we have used to compute the discon- 
tinuities. Let us just try outline the method for a one- dimensional integral such as Eq. (4.21). 
Basically, there are exactly three phenomena to keep track of: 

• A pole of the integrand makes a loop around an integration endpoint. 

• A branch cut endpoint of the integrand makes a loop around an integration endpoint. 

• The value of the integrand at an endpoint undergoes monodromy. 

By keeping track of these phenomena, we can calculate the symbol of an arbitrary integral 
of the form Eq. (4.21). The procedure is systematic and entirely algebraic, and in principle 
could be automated. A simple yet nontrivial example illustrating it is given in Appendix B. 

We have applied this algorithm to the transcendentality-two integrand in Eq. (4.21). At 
the worse point the procedure generated no more than order 50 terms for a given j. Still, 
several significant simplifications occured, and the resulting symbol is reproduced in Appendix 
A. For the reader's convenience, it is also attached with this submission to the arXiv in the 
form of a Mathematica notebook. Nicely, by adding suitable telescopic terms in j, we found 
that the symbol can be expressed as a sum over j of a general term A, which is both integrable 
manifests the correct symmetry under parity combined with the reflection 2 o i. 

The parity-odd part of A turns out to match the 3-mass "easy" six-dimensional hexagon 
computed recently in [52-55] (this can also be seen directly from Eq. (4.21)). It would be 
nice to find a similar interpretation for the parity even part. 
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Let us summarize. The differential of the remainder function is 

dR^ = Cijdlog{i-lii+lj) (4.26) 

where 

a,, = cQ + cQ + cf^ + C^- (4.27) 

The functions C^^) and C(^) are given analytically in Eqs.(4.9) and (4.25). The function C^^-* 
admits the integral representation (4.21) and C*^^-* is obtained from it by symmetry. We did 
not attempt to express C^^-* explicitly in terms of polylogarithms, but its symbol is given in 
Appendix A. 

5 Checks on the two-loop result 

In this section we discuss some consistency checks which, in our opinion, establish beyond 
reasonable doubt that the obtained expression for the two-loop remainder function is correct. 

5.1 Hexagon case 

We have verified that the symbol given in Appendix A reproduces that of the compact analytic 
formula given in [1] when n = 6. 

5.2 Integrability 

The differential in Eq. (4.26) is only meaningful if it is consistent with = 0. This is what 
we call the integrability condition. We did not attempt to verify this at the level of the 
integral representation, but we did verify, using a computer, that it was true at the level of 
the symbol. This establishes, at least, that the maximal functional transcendentality part of 
our differential is the differential of a function. This property holds for the whole C, but not 
for the individual C^"^^ It would be nice to understand how integrability derives from the 
supersymmetry algebra in general. 

5.3 Physical discontinuities 

The discontinuities of a scattering amplitude have a clear physical interpretation. For in- 
stance, the leading branch cuts should be unitarity cuts, and the location of a branch point 
should be some physical threshold. In a theory with only massless particles, these thresholds 
occur when a separation [xi — Xj)'^ becomes null. Therefore, the first entry of the symbol 
is expected to be always of the form (ii+ljj+l). We can keep going. The branch cuts 
themselves have branch cuts. While their physical interpretation is more delicate, the kine- 
matical considerations which dictate the location of the endpoints can involve at most one 
extra channel and therefore should be similar from one- to two- loops. The endpoints of higher 
discontinuities, however, may in principle depend on three or more channels at a time, and 
could be very complicated conditions. 
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Thus, specifically the first two entries of the symbol are expected to contain only patterns 
of the form which appear in a general 1-loop amplitude, like Li2(l — u), or involve two inde- 
pendent cuts like \og{ii+ljj+l) log{kk+lmm+l) . (These expectations were also formulated 
in [15].) We believe that a general proof of this statement is possible by considering Feynman 
integrals [51]. Pleasingly, we find that our symbol fulfills precisely these expectations. 

5.4 Collinear limits 

The remainder function i?„ is designed to have simple collinear limits [5] . Let us consider the 
simplest such limit, where two adjacent edges of the polygon become parallel. In this limit 
the n-gon is reducing to a (n— l)-gon, losing three degrees of freedom. The approach to the 
limit is therefore characterized by three parameters. We find convenient to use the following 
parametrization of momentum twistors: 

Zn = Zn-l + e{Zn-2 + TZi) + e'Z2 (5.1) 

where e' ^ e. The parameters e and e' regulate the approach to the collinear limit, while r 
describe the relative length of the collinear segments n and n—1. For the remainder function, 
the dependence on all three parameters must drop out in the limit: 

lim RniZi, . . . ,Zn) = R'n-AZi, . . . ,Zn-l). (5.2) 

This imposes nontrivial constraints on the differential of R, namely 

lim dR = dR' (5.3) 

e,e'-S>0 

which we now wish to verify. Our expression for the differential is a sum of terms of the form 
Cijdlog{i—lii+lj). The log derivative of the small parameter e' multiplies {Cn-i,i + Cn,n-2)] 
this combination must vanish. Similarly, from canceling the log derivative of e and using the 
exact sum rule Eq. (3.10), we conclude that {Ci^n-i + Cn-2,n) vanishes in the limit. From 
the pole at dr/r at r = 0, which multiplies {Cn-2,n — Cn-i,i), we conclude that both must 
vanish separately. Similarly from the pole at r = oo. (Strictly speaking we can only conclude 
that these quantities vanish separately when r = or oo, but this provision does not seem to 
be realized for the 2-loop amplitude.) There remains the finite parts to match on both sides, 
which produce a large number of relations which can be summarized very simply 

limCjj = C-j, i,j 7^ n— 1,71, (5.4a) 
limCi,n-i + Ci^n = C'in-i, i / n-l,n, (5.4b) 
limCn-i,j + Cn,j = C'n-i,ji J / n-l,n, (5.4c) 

where the C pertain to the (n— l)-gon. The extra conditions, lim C„_i^i = and lim Cn,n-2 = 
0, can be viewed as special case of these. 

As above, we have verified using a symbolic manipulation package that the limits (5.4) 
work out at the level of the symbol (up to n = 10, which is generic). 
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5.5 Two-dimensional kinematics 

As another important check on our results, we have compared with existing results in two- 
dimensional kinematics, which were given for an arbitrary number of points in [3, 4]. In this 
case we were able to fully check the integral representation, Eq. (4.21), not only its symbol. 

In two-dimensional kinematics, the number of particles is taken to be even and particles 
with odd and even labels are distinguished, 

Z2._i ^ 0, 0) , Z2i ^ (0, 0, . (5.5) 

Four-brackets with two odd and two even labels factorize into products of two-brackets, and 
all other four-brackets vanish. Provided that Ceven.odd 0, Eq. (3.15) becomes 

dRn^ {Cij + Cj,i)dlog{ij) + Yl iCi,j+Cj,i)dlog{ij). (5.6) 

i<j i<j 
both odd both even 

We would like to evaluate the functions Cij in these kinematics. First of all, the vanishing 
of Ceven,odd is easy to see. For odd i, we have that: 

(2i-l)(2i+l)(H)(3i) 



form which C^^^ and C^^^ trivially vanish. In two-dimensional kinematics, we write X = Z2AA 



(2') 

where A = Z\ — tZ^^ in Eq. (4.21). It is then also easy to see that C2 / vanishes for i odd. 



because the measure becomes dx log — )• dA log |^ = 0. 

Thus C2,j is only nonzero for even i. These are slightly messy to compute, because the 
way we have organized our computation is not really tailored for the special simplifications 
that occur in the special kinematics. We give some details of this computation in Appendix 
C. In the case of the octagon, the final result of that Appendix is simply 

C2,4 = -logni,3,5,7log-U3,5_7,i logU4,6,8,2, 
C2,6 = logni,3,5,7logn3,5j,i log(n2,4,6,8^i4,6,8,2) 

plus cyclic. Equation (5.6) can then be easily integrated: 

= -2 log ui,3,5,7 log ■"3,5,7,1 log ^2,4,6,8 log ^4,6,8,2 + Constant. (5.8) 
10 



This agrees precisely with [4]. The constant cannot be determined by our methods; its 
9 



value. Constant = — ^f-, may be fixed from triple collinear limits and the hexagon [4].^^ 



11 



We are factoring out a 2-loop prefactor where = ^let^ " ' ^'^ '^^^ times that of [4 

The author thanks Claude Duhr and Vittorio Del Duca for explanations regarding this point. 
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6 Conclusions 



We have considered, in planar = 4 super yang-Mills, a generalization of scattering ampli- 
tudes and conjectured that it obeys simple Ward identities under dual superconformal symme- 
tries. These Ward identities follow from naive field theory manipulations, and the non-obvious 
part of the conjecture is that they should not be spoiled by any sort of divergences, which 
arose previous treatments of superconformal symmetry on scattering amplitudes. Somewhat 
analogously to the Ward identities employed in the classic work of Parke and Taylor [57] 
on MHV amplitudes nearly three decades ago, this allows a difficult computation involving 
gluons to be replaced with a less difficult ones involving fermions. In the present case, the 
computation with fermions is related to a derivative of the computation with bosons. 

Exploiting these Ward identities we obtained an integral representation for the total 
differential of the two-loop general n-point MHV remainder function. The result agrees with 
all previously known analytic two-loop results, and extends them. Its symbol is reproduced in 
Appendix A. We have really only used N = 1 supersymmetry, and our results are thus valid 
for null polygonal Wilson loops in any N = 1 superconformal theory. The next natural step 
would be to integrate this expression. It would also be worthwhile to compare it numerically 
with the results [58]. With some work, we believe the method could be realistically expected 
to generate new results for three-loop MHV or two-loop NMHV amplitudes [59]. 

The generalization is based on the introduction of new variables which render exact all 
the dual supersymmetries, e.g. those which act naturally on Wilson loops. This includes 
a chiral half of the spacetime superconformal symmetries. An important question, which 
remains to be addressed, is whether the other chiral half also becomes exact. This would 
imply exact Yangian symmetry for the same price. This would also suggest a generalization 
for scattering amplitudes in superconformal theories other than N = A. 

The methods we have used should commute with the Wilson loop OPE [13], and with 
structural constrains derived from it [60]. It might be fruitful to combine the two approaches. 

One may ask whether general constraints can be derived regarding the structure of am- 
plitudes, coming from the existence of this new object. Our results provide evidence that 
this is the case. Of course, in non-chiral superspace any function can be supersymmetrized, 
so supersymmetry alone does not provide much of any constraint. However, we have found 
empirically that individual 6, 9 components of the amplitude at two-loop have unexpectedly 
simple leading singularities, all equal to "1" in the natural normalization, which in turns 
implies a very simple form for the last entry of the symbol at the same loop order. This 
property begs for an explanation, which we hope will help unlock the structure of the symbol 
at higher loops. 
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A Symbol of Cij 

The differential of the n-point function is expressed as 

dRn = Cijd\og{i-ln+lj) 

where C2,i is the sum of the four contributions 

j— 1 n+l 



(A.1) 



a 



(1) 



2,i 



iog«2,i-i,i,i X 

j=2 k=i 

i-2 

<^2? = XA(l,2;j-l,j;i-l,i) 

j=4 



Liofl - u 



■j,k,k-l,j+l. 



+ log 



X 



log 



^lk-1 



Cf} = J2 A(2,l;j,i-l;i,i-l), 



a 



(4) 



j=i+2 

-2Li3(l - -) - Li2(l - -) log« - J log3 u + ^ logn, 
n n 6 6 



(A.2) 



and other Cij arc obtained by cyclic symmetry. In the first line, Xj^i = X2 when j = 
and Xk-i = xi when k = i, and in the last line, u = U2,i-i,i,i- The symbol of A is 

SA{l,2;j-l,j-i-l,i) 

(u+l(2)n(j))(23ii) 



= S[hir,l,2;j-l,j) 



{j-ljj+U){123j){23ii+l) 



{{ii+l) (i-li)) 



+ 



^<;n\n(^ • A-^\n(^ -n,-.. A] ( (123^)0^-lii+12)(23^j) \^ 0i+l(2)n(i))(u+lii+l) 

20[l^l2;^i %2,l,i-lj 1^121-1- %2,l,jjJ 09 \^(i23j>0-ljj+H){23M+l)^ {2ijj+l>(13(2i-H)n(2jj+l)) 

+i<S[Li2(l - nj,j_i,i,2) - Li2(l - (g) (^^ 



(12i-H)(23ij) 



{123i>{i-lii+li>{23i-li> 

-^ 1 9rT ion - 7/0 ■ 1 ■ , ^1 (Si (iJ+i(2)n(i))(i-H+i(»23)nfaj+i)) 

+ 2>^[l^l2li - U2,t-l,t,l}\ (2iii+l>(13(2i-li)n(2ii+l)> 

-^l'^^np■7/•• , •nlnp-7/-o, • , 1 fc, (23^) (ij+l(2)n(i))(13(2t-H)n(2ii+l)) 
+ 2'^LiOg«J,«-l,«,2iOg«j,2,l,i-lJ OS l^(123i> ) (2iji+l>(23i-li>(i-li+l(j23)n(iji+l)> 

V -((ii+l) ^ (j-ii)) 

(12ij)(23i-li) 



(ij+i(2)n(i))(i-ii+i(i23)nfaj+i)) 

{2iii+l>(12ij+l> 



+ 5[/5(l;i-l,i;j-l,j) 



+ 5[log log ^X2,i-l,i,l 



(12z-H)(23zj) 

(i-ii+i(ii2) n (in+i))(i23i)(23i-ii) 



(A.3) 



(123i)0-ljj+H)(12z-H)(23zi+l) 

The factors of ^ cancel telescopically in the sum over j, and there are no ^ in front of anything 
in the full symbol of the amplitude (e.g., inside the big parenthesis, only the squared factors 
do not telescope away). The symbol could be written more succintly by exploiting these 
telescopic cancellations; this particular presentation makes the individual term A integrable 
and parity covariant. is the "pentagon integral" 

/5(X; l,2;i-l,z) = Li2(l - ttx,i,2,i) -Li2(l -ux,i,2,i-i) + Li2(l - «2,i-i,i,i) 

+Li2(l - ux,i,i-i,i) - Li2(l - ux,i,i-i,2) + logMx,i,2,ilogttx,i,i-i,i- (A.4) 
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B Procedure for computing the symbol 



In this section we wish to illustrate the method for computing the symbol outlined in subsec- 
tion 4.5, on a simple example; we hope to expand more on the method elsewhere. 

Example = / . ; log(l + x)=Li2(l-b)-Li2(l-a). (B.l) 

Jq {x + a){x + b) 

We want to compute its symbol without actually doing the integral. The first step in the 
method is to find a function dual to the boundaries of the integration region: /(x) = x/d, 
where d is an arbitrary nonzero constant. This function is characterized by having a pole on 
the endpoint at infinity and a zero at the other endpoint. The three phenomena referred to 
in the main text then contribute as follow. The first phenomenon always contribute, for a 
pole located at point —a, an amount /(—a) tensored with the residue of the pole: 

^®(l-a). (B.2) 

(Signs can be dropped for entries inside the symbol.) The pole at x = b contributes 

_ b 
~ d 

The second phenomenon contributes /(—I) tensored with an integration which ends at the 
branch point 1: 

The third phenomenon replicates an existing endpoint, and moves the leftmost entry of the 
symbol of the integrand, evaluated at the endpoint, outside the integral: 

Jo{x + a){x + b) 

Note that in this case the first entry is 1, so this contribution vanishes. There is also the 
endpoint at infinity. In the present example, phenomena 2 and 3 occur simultaneously at 
that point. A simple and systematic way to untangle this, is to separate the two endpoints 
by chopping the integral at A, so that / = xA/d{x — A), and then take A — t- 00: 

A ^ (ixia -b) A ^ /-"^ dx{a -b) , ^ dx{a - b) 
{l-A)(g)S / ^— -IX-Z,®'^ / 7 -r^^d(g)S / (B.6) 



(B.3) 



{x + a){x + b) d J {x + a){x + b) J {x + a){x + b) 

The integrals we have written down have only one endpoint and, correspondingly, would seem 
to not make sense. A nontrivial fact, which will always be guaranteed, is that the endpoints 
organize into pairs. Here the contributions of the second and third type add up to 

, ^ dxia-b) , 1-a ,^ , 

J I {x + a)[x + b) 1-0 
A second nontrivial fact, which can be used either as a check on the computation, or as a 
way to remove all but (B.2) and (B.3), is that the dependence on d cancels in the total: 

^[/example] = a ® (1 - a) - 6 ® (1 - 6) (B.8) 

in agreement with the explicit result of doing the integral. 
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C Two-dimensional kinematics 



In this section we complete the computation of the function Cij in the kinematics (5.5). The 
function l5{X; that enters Eq. (4.21) vanishes when both i and j are even (this 

can be seen directly, pre-integration, from the vanishing of the numerator in Eq. (4.6)). For 
j odd it is pure logarithms. (Passing from forms with dilogarithms to forms with logarithms 
is easy, if one uses the symbol, and fixes vr^/G ambiguities by using special points.) After 
collecting some telescopic terms in Eq. (4.21) we find "only" 

f' , (-4i-l) (Az-l)(i+lj) (2j+l)(j-H) 



{Aj) ^ {Ai+l){i-lj) -(2j-l)(j+H) 



j=5 odd 

(^^-1), {Ai-l){jl)^ (2j+l){j-li) 



j=5 odd 

r3 



{Ai-l) / . .(lA)(i+13), ^. ,(l^)(i-13), ^. ,{lA){i+li-l) 



+ A WTT) r^^z+l^)(13)) - ^^^((^-l^)(13)) - (.+1^)(H-1) 

+ log!f^'!!^"''!log!'^'^w^.^''n - /^.log^Li.(-^^^>^^-^^+^>^ 



(Ai-l)(i+ll) ^ (i+13)(li-l)y Ji ^ {A3) ^'{Ai-l){3i+l)'' 

This looks a mess but it is not so bad. The first two lines sum up to simple logarithms. The 
remaining dilogarithms are a bit vexing, but they can be removed by a simple integration by 
parts; this sort of simplification is specific to the special kinematics. Therefore 

(2) _ {3i-l){i+lj) (3i-l)(jl) (2j+l)(j-H) 



Co / = > log -, n r log -, — n r log 



,=5„dd (3i+l)(^-lj) (3i)(^-ll) "(2j-l)(j+H) 

+ diog — ^log log ^ - 2 log^d + j • (c.i) 

The second line is still a bit ugly. However, beautiful things happen when we add to it C(3\ 
which contains the same term, and C^^^: the sum is simply loguLi2(n). This cancels a term 
from C^^\ which we can then write compactly 



62^] + log 'uLi2('u) = log u X 



/ E^=5odd ELi+2even log -"S.fc+l.fc-lj logUij^ij+i^k 

+ E}=4even Efc^i+Sodd ^OgU2,k+l,k-l,j log Uj+i j_i | . (C.2) 

y + log(l - u) logU2,2,i,i 



The cross-ratio U2,2,i,i is not a typo, it means we are somewhat abusing notation to produce a 
more compact expression: the ill-define brackets (22) cancels out e.g. between the second line 
with k = n+1 and the third line. Our final result for C2,j is the first line of Eq. (C.l), plus an 
equivalent term for C^^\ plus Eq. (C.2). For the octagon this is discussed in the main text. 
We have verified symbolically that twice this result agrees with the derivative with respect 
to (2z) of the general formula in [3] (up to n = 20, which is amply generic), confirming the 
results of section 4 in the case of special kinematics. 
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